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Abstract 

Two-body charmless hadronic B decays involving a tensor meson in the final state are studied 
within the framework of QCD factorization (QCDF). Due to the G-parity of the tensor meson, 
both the chiral-even and chiral-odd two-parton light-cone distribution amplitudes of the tensor 
meson are antisymmetric under the interchange of momentum fractions of the quark and anti- 
£N) ■ quark in the SU(3) limit. Our main results are: (i) In the naive factorization approach, the 



decays such as B — > -k~ and B — > K^tt with a tensor meson emitted are prohibited 
owing to the fact that a tensor meson cannot be created from the local V — A or tensor current. 
Nevertheless, they receive nonfactorizable contributions in QCDF from vertex, penguin and hard 
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spectator corrections. The experimental observation of B~ — > K^tt - indicates the importance of 
nonfactorizable effects, (ii) For penguin-dominated B — > TP and TV decays, the predicted rates 
in naive factorization are usually too small by one to two orders of magnitude. In QCDF, they are 
enhanced by power corrections from penguin annihilation and nonfactorizable contributions, (iii) 
The dominant penguin contributions to B — > K^if* ) arise from the processes: (a) b — > sss — > st] s 
and (b) b — > sqq -4 qK% with r\ q = (uu + dd) / \/2 and rj s = ss. The interference, constructive 
for K^rj' and destructive for K^r], explains why T(B — > K^rf) » T(B — > K^rj). (iv) We use the 
measured rates of B — > K^ito, 4>) to extract the penguin-annihilation parameters p T j^ and p]^ T and 
the observed longitudinal polarization fractions fL(K%ui) and f^K^) to fix the phases and 
■ ( v ) The experimental observation that /t//l <C 1 for B — > if|(1430)(/), whereas /t//l ~ 1 
for B — > ET|(1430)o; with fx being the transverse polarization fraction, can be accommodated 
in QCDF, but cannot be dynamically explained at first place. For penguin-dominated B — > TV 
decays, we find /^(K^p) ~ /l^K^uj) ~ 0.65, whereas fi{K*f2) ~ 0.93. It will be of great 
interest to measure /l for these modes to test QCDF. Theoretically, transverse polarization is 
expected to be small in tree-dominated B — > TV decays except for the a^p , a 2 P + 1 K^K*~ 
and K^K* modes, (vi) For tree-dominated decays, their rates are usually very small except for 
the 02(71"^, «2 (7r - ,/9~ ) and /2(vr _ ,p _ ) modes with branching fractions of order 10~ 6 or even 
bigger. 
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I. INTRODUCTION 



In the past few years, BaBar and Belle [l- 14] have begun to measure several charmless B decay 
modes involving a light tensor meson T in the final states with the results summarized in Table 
HI From the theoretical point of view, the hadronic decays B — > TM with M = P,V,A are of 
great interest for two reasons: rate deficit and polarization puzzles. First, these decays have been 



studied in the naive factorization approach 1514231] . The predicted rates are in general too small by 



one to two orders of magnitude. This implies the importance of l/rrif, power corrections. Since the 
nonfactorizable amplitudes such as vertex and penguin corrections, spectator interactions cannot 
be tackled in naive factorization, it is necessary to go beyond the naive factorization framework. 
The theoretical frameworks suitable for this purpose include QCD factorization [3], perturbative 



QCD (pQCD) [25(] and soft-collinear effective theory (SCET) [2J 

Second, it is known that an unexpectedly large fraction of transverse polarization has been ob- 
served in the penguin-dominated B —> VV channels, such as B ^ (pK* , pK* , contrary to the naive 



expectation of the longitudinal polarization dominance (for a review, see 27]). However, while the 
polarization measurement in B — > wX|(1430) indicates a large fraction of transverse polarization 
fx (see Tabled]), the measurement in B — )■ 0X|(143O) is consistent with the longitudinal polariza- 
tion dominance. Therefore, it is important to understand why /t//l *C 1 for B — > 0X|(143O), 
whereas fr/fz ~ 1 for B — )■ wAT|(1430), even though both are penguin-dominated. The polariza- 
tion studies for B — > TV, TA, TT will further shed light on the underlying helicity structure of the 
decay mechanism. 

In the present work we shall study charmless B — > TM decays within the framework of QCD 
factorization. One unique feature of the tensor meson is that it cannot be created from the V — A 
or tensor current. Hence, the decay with a tensor meson emitted, for example, B~ — > X|°7r _ , is 
prohibited in naive factorization. The experimental observation of this penguin-dominated mode 
with a sizable rate implies the importance of nonfactorizable effects which will be addressed in 
QCDF. 

The layout of this work is as follows. We study the physical properties of tensor mesons such as 
decay constants, form factors, light-cone distribution amplitudes and helicity projection operators 
in Sec. 2 and specify various input parameters. Then we work out in details the next-to-leading 
order (NLO) corrections to B — > TP, TV decays in Sec. 3 and present numerical results and 
discussions in Sec. 4. Conclusions are given in Sec. 5. Appendix A is devoted to a recapitulation 
of the ISGW model. Decay amplitudes and explicit expressions of helicity-dependent annihilation 
amplitudes are shown in Appendices B and C, respectively. A mini review of of the rj — rf mixing 
is given in Appendix D. 

II. PHYSICAL PROPERTIES OF TENSOR MESONS 
A. Tensor mesons 

The observed J p = 2+ tensor mesons / 2 (1270), $(1525), a 2 (1320) and AT 2 *(1430) form an SU(3) 
1 3 i-2 nonet. The qq content for isodoublet and isovector tensor resonances is obvious. Just as the 
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TABLE I: Experimental branching fractions (in units of 10 6 ) and the longitudinal polarization 
fractions fi for B decays to final states containing a tensor meson. Data are taken from 



Mode 




B 


/x 


Mode 




B 


h 


B(B+ - 


■»■ iq(1430)+w) 


21.5 ±4.3 


0.56 ±0.11 


B(B° - 




10.1 ±2.3 


0.45 ±0.12 


B{B+ - 


* K$(1430)+</>) 


8.4 ±2.1 


0.80 ±0.10 


B{B° - 


+ if 2 *(1430)V) 


7.5 ± 1.0 


n qm+0.059 
u.jui_ 069 


B{B+ - 


+ K 2 *(1430)+ ?? ) 


9.1 ±3.0 




B(B° - 


> #2* (1430)%) 


9.6 ±2.1 




B{B + - 


* K* 2 (U30)+ V ') 


28.0±g;B 




B(B° - 


> ir|(1430) r/') 


13.7±|;f 




B(B+ - 


* K 2 *(1430)°vr+) 


5.6±?;| 




B(B° - 


* iq(1430)+7r-) 


< 6.3 




B(B+ - 


> K 2 *(1430) X+) 


< 1.1 




B{B° - 


* ir|(1430) vr ) 


< 4.0 




B(B+ - 


* $(1270)K+) 


1 Ofi +a28 

l.UD_ 29 




B(B° - 


+ / 2 (1270)tf°) 


O 7+1-3 
z - ' -1.2 




B(B + - 


+ / 2 (1270)vr+) 


J-- ' -0.49 












B(B+ - 


+ $(1525)K+) 


< 7.7 a 












B(B + - 


■»■ a 2 (1320)°K + ) 


< 45 













a From the measurement of B(B+ -> $(1525)°A'+ -> K+K+R-) < 3.4 x 1CT 6 [6j. 



rj-rj' mixing in the pseudoscalar case, the isoscalar tensor states /2(1270) and $(1525) also have a 
mixing, and their wave functions are defined by 



/a (1270) 
$(1525) 



1 

i 

71 



(/ 2 " + $ d )cos0 /2 + $sin0 /2 
(# + / 2 d )erin0 A -/|coe0 /a 



(1) 



with f% = uu an d likewise for f!j' s '■ Since tttt is the dominant decay mode of /2(1270) whereas 
$(1525) decays predominantly into KK (see Ref. [28l ) , it is obvious that this mixing angle should 
be small. More precisely, it is found that 0/ 2 = 7.8° 29] and (9 ± 1)° 0. Therefore, $(1270) is 
primarily an (uu + dd)/y/2 state, while $(1525) is dominantly ss. 

For a tensor meson, the polarization tensors ef^ with helicity A can be constructed in terms of 



(2) 

(3) 
(4) 
(5) 



the polarization vectors of a massive vector state moving along the z-axis 
e (0)*" = (P 3) 0, 0, E)/m T , e(±l)*" = (0, T l, +i, 0)/V2, 



and are given by 



£ (±2) 
£ (±1) 
£ (0) 



e(±l)"e(±l)" 

T 

2 
T 



e(±l)^e(0) 1 ' + e(0fe(±l)H, 



e +1 



■ + ^(0)^(0)". 



The polarization e^, can be decomposed in the frame formed by the two light-like vectors, z^ 
and p v = P v — z v m^/(2pz) with P v and tut being the momentum and mass of the tensor meson, 



respectively, and their orthogonal plane 3l|, |32|. The transverse component that we use thus reads 



e (A) z v 



c fiv 41 



(A) u 



c jj,u *• 



,(A) a „v 
pz 



in 



T 



2pz 



(6) 
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The polarization tensor satisfies the relations 



0, 



p f^ v 



The completeness relation reads 

E e< $ { e( P X J)* = \ M W M VU + l -M pa M vp - X -M pv M prj , 

A 

where M pv = g pv - P p P u /m^. 



(7) 



(8) 



B. Decay constants 

Decay constants of the vector meson are defined as 



(V(P,e)\q 2 ^ gi \0) 
(V(P,e)\q 2 a^ qi \0) 



(9) 



Contrary to the vector meson case, a 3 P 2 tensor meson with J PC = 2 ++ cannot be produced 
through the local V — A and tensor currents. To see this, we notice that 



(T(P,\)\V,\0) = ae< x >P v + be< x >P li = 0, 
(T(P,X)\A^\Q) = e pupf7 P u e^* = 0, 



(10) 
(11) 



where use of Eq. ([7]) has been made. Nevertheless, a tensor meson can be created from these local 
currents involving covariant derivatives: 



<T(P, A)|<V(0)|0> = f T m 2 T e;i x \ 



(T(P,A)|J^(0)|0) 



timT(e$*P u -eM*P f 



(12) 



where 



<V(°) = 2(« 1 (% i ^ ( '2(0) + 9i(0)vB /1 ?2(0) 



J ia(°) = 9i(0)<V«A*g2(0), 



(13) 



and D P =D^ - D p with D p =d, +ig s A a p X a /2 and D^=d» -ig s A<l l X a /2. 

The decay constant fr of the tensor meson has been estimated using QCD sum rules for the 
tensor mesons /2(1270) 33J and X|(1430) 34] and the tensor-meson-dominance hypothesis for 



/2(1270) [33j, |35|, [36| . The previous sum rule predictions are [33j, [34| 1 

// 2 (1270)G" = 1 GeV) ~ 0.08m /2(1270) = 102 MeV , 
/jk(143O)0* = 1 GeV) ~ (0.10 ± 0.01) m K * (1430) = (143 ± 14) MeV . 



(14) 



Recently, we have derived a sum rule for /t(a*)/t (p) and revisited the sum-rule analysis for /t(aO- 
Our results of /t and /jr for various tensor mesons are shown in Table HVl below 40]. Our sum 



1 The dimensionless decay constant fx denned in 33|, |34[ differs from ours by a factor of 2ttit- The factor 
of 2 comes from a different definition of D p there. 
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rule results are in good agreement with [33f] for // 2 (i27o)> but smaller than that of [34( for fK*(i±3Q)- 
The decay constants for /2(1270) and /2(1525) also can be extracted based on the hypothesis of 
tensor meson dominance together with the data of r(/2 —> tttt) and IX/2 —> KK). We found that 



4C] 

// 2 (1270) - (0-085 ± 0.001)m /2(1270) = (108 ± 1) MeV , 

// 2 (1525) - (0-089 ± 0.<m)m% (1S2S) = (136 ± 5) MeV . (15) 
They are in accordance with the sum rule predictions shown in Table lIVl 



Form factors 



Form factors for B —> P,V,T transitions are defined by 42h44 ] 



(P(P)\V,\B(p B )) 
(V(P,X)\V,\B( PB )) 



2 



m T3 — m. 



% ) if P (S 2 ) + 



m r — mi 



{V{P,\)\A„\B{j, B )) = 2m v 



m B + my 



s^e^p a B P^V BV (q 2 



Y-^-q^A^ \q 2 ) + (m B + my) 



PB 



(T(P,X)\VJB(p B )) 



<r(P,A)|^|S(pB)) = 2m T 



m B + my 
2 

-1 

m B + tut 

e (A> . 



m n — m 



P/j, + (pb) 



-^e^ )P %P^V BT (q% 



v 



A BV {q^ 



j^q^iq 2 ) + (m fl + m T ) 



m B + m-r 



^ + (Pb) 



m B — 777,5 



^ V (<Z 



A BT {q% 



A BT ^ 



(16) 



where = (p B — P)^ and er-X = ^(x^P Bv l m B- Throughout the paper we have adopted the 
convention e 0123 = — 1. 

In the Isgur-Scora-Grinstein-Wise (ISGW) model [4^], the general expression for the B — > T 
transition is parametrized as 

(T(P,X)\(V-A) fl \B(p B )) = iKq 2 )e^ pa e* va p Ba {p B + Pyq° - k{q 2 )^ v p B 

- b + {q 2 y ap p%p B {p B + P), - b-(q 2 y a pp a B p B q„ (17) 

where the form factor k is dimensionless, and the canonical dimension of h, 6+ and 6_ is —2. The 
relations between these two different sets of form factors are 



V BT (q 2 ) = m B {m B + m T )h(q 2 ), A B1 (q 



BT/2\ 



m B 



m B + tut 



k(q 2 



(18) 



A BT (q 2 ) = -m B (m B + m T )b + (q z ), A* 1 (q 



BT I „ 2\ m B \,2/J2 



2rriT 



[k 2 {q') + (m B - m' T )b + {q') + q'b.(q% 



The B — > T transition form factors have been evaluated in the ISGW model 



45J and its improved 



version, ISGW2 [46|], the covariant light-front quark model (CLFQ) |47j], the light-cone sum rule 
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(LCSR) approach [48J], the large energy effective theory (LEET) [49H51I] and the pQCD approach 
44|. In LEET, form factors are evaluated at large recoil and all the form factors in the LEET limit 



to be s pec ified below can be parametrized in terms of two independent universal form factors (j 
and C|| B : 



V BT (q 2 ) 
A BT {q 2 ) 



\'Pi 



Wt\ 
m T 

Wt\ 



m B J 



rriT 



2Ejn 



rriB + rriT 
1+ mr\ 

\P T \ \ m B / 

where Et is the energy of the tensor meson 



atf ) - ^tf> 



In the LEET limit, 



E T ,m B > m T , Aqcd- 



(19) 



(20) 



(21) 



Using the recent analysis of tensor meson distribution amplitudes [40J], one of us (KCY) has 
calculated the form factors of B decays into tensor mesons using the LCSR approach 48]. The 
LCSR results are close to LEET and pQCD calculations. 

The B — > a2(1320), f2 q = {uu + dd) /v2, 7^1(1430) transition form factors calculated in various 
models at the maximal recoil q 2 = are summarized in Table HH The ISGW model [2] is based 
on the non-relativistic constituent quark picture. In general, the form factors evaluated in the 
ISGW model are reliable only at q 2 = = (tub — wit) 2 \ the maximum momentum transfer. The 
reason is that the form-factor q 2 dependence in the ISGW model is proportional to exp[—(q 2 n — q 2 )] 
and hence the form factor decreases exponentially as a function of (q^ — q 2 ) (see Appendix A for 



details). This has been improved in the ISGW2 model 46|] in which the form factor has a more 



ravior at large {q 2 ^ — q 2 ) which is expressed in terms of a certain polynomial term. As 



201 ] . form factors are increased in the ISGW2 model so that the branching fractions of 



realistic be : 
noticed in 

B — > TM decays are enhanced by about an order of magnitude compared to the estimates based 
on the ISGW model. 

The CLFQ model is a relativistic quark model in which a consistent and fully relativistic treat- 
ment of quark spins and the center-of-mass motion is carried out. This model is very suitable 
to study hadronic form factors. Especially, as the recoil momentum increases (corresponding to 
a decreasing q 2 ), we need to start considering relativistic effects seriously. In particular, at the 
maximum recoil point q 2 = where the final-state meson could be highly relativistic, it is expected 
that the corrections to non-relativistic quark model will be sizable in this case. 

The CLFQ and ISGW2 model predictions for B — > T transition form factors differ mainly in 
two aspects: (i) when q 2 increases, h(q 2 ), 6+(<? 2 ) and b-(q 2 ) increases more rapidly in the former 
and (ii) the form factor k obtained in both models is quite different, for example, k BK 2 (0) = 0.015 
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TABLE II: B -»• T transition form factors at = evaluated in the ISGW2, CLFQ, LCSR, 
LEET and pQCD models. The CLFQ results are obtained by first calculating the form factors 
h(q 2 ),b+(q 2 ) and b^(q 2 ) using the covariant light-front approach and k(q 2 ) from the heavy quark 
symmetry relation Eq. (|22p and then converted them into the form- factor set V(q 2 ) and ^4o,i,2(<? 2 )- 
To compute the form factors in LEET we have applied C±(0) = 0.28 ±0.04 and C||(0) = 0.22 ±0.03 . 
LCSR and pQCD results are taken from [48] and [3], respectively. 



ISGW2 CLFQ 



LCSR 



LEET 



pQCD 



n 1s +o.()5 

11+ 003 
u,11 -0.03 

0.081°;°* 



ISGW2 CLFQ 



LCSR LEET 



A o 

A.Ba 2 
A 2 

A Bh " 

A 
A 2 
A 



pQCD 



Timr 

-0.04 
+0.02 
-0.01 

0.25 0.20 ±0.04 0.13 ±0.02 0.13 jl ,;^ 
0.19 0.10 ±0.02 0.13 ±0.02 0.04+°;^ 
0.23 0.25 ±0.04 0.15 ±0.02 O.lSj:^ 
0.21 0.05 ±0.02 0.14 ±0.02 0.08^;^ 



yBa 2 
A Ba 2 

yBj 2q 

A Bf 2q 
IfBKZ 



,dk: 



0.32 0.28 0.18 ±0.02 0.18 ±0.03 

0.16 0.21 0.14 ±0.02 0.13 ±0.02 

0.32 0.28 0.18 ±0.02 0.18 ±0.02 

0.16 0.21 0.14 ±0.02 0.12 ±0.02 

0.38 0.29 0.16 ±0.02 0.21 ± 0.03 

0.24 0.22 0.14 ±0.02 0.14 ±0.02 



0.21 



+0.06 
-0.05 



13+ 004 



0.20 
0.14 
0.20 
0.14 
0.27 
0.22 



0.24 0.21 ±0.04 0.14 ±0.02 0.18 
0.19 0.09 ±0.02 0.13 ±0.02 0.06 



TABLE III: B — )■ T transition form factors obtained in the covariant light-front model and fitted 
to the 3-parameter form Eq. (|23p . 



F 


F(0) 


a 


b 


F 


F(0) 


a 


b 


yBa2 


0.28 


2.19 


2.22 


A Ba 2 


0.24 


1.28 


0.84 


Af a2 


0.21 


1.38 


0.47 


jy3a 2 


0.19 


1.93 


1.69 




0.28 


2.19 


2.22 




0.25 


1.37 


0.95 


A B hq 


0.21 


1.39 


0.46 




0.19 


1.93 


1.69 


V BK 2 


0.29 


2.17 


2.22 


.BK* 

A 2 


0.23 


1.23 


0.74 


A* K * 


0.22 


1.42 


0.50 


.BK* 

A 2 2 


0.21 


1.96 


1.79 



in the former and 0.293 in the latter. Indeed, it has been noticed [47j that among the four B — >• T 
transition form factors, the one k(q 2 ) is particularly sensitive to /3t, a parameter describing the 
tensor-meson wave function, and that k(q 2 ) at zero recoil shows a large deviation from the heavy 
quark symmetry relation. It is not clear to us if the very complicated analytic expression for k(q 2 ) 
in Eq. (3.29) of [47] is complete. To overcome this difficulty, it was pointed out in 47] that one 
may apply the heavy quark symmetry relation to obtain k(q 2 ) for B — >■ T transition 



k(q 2 



mBrriT 1 + 



h{q 2 )- l -b + (q 2 )+ l -b^q 2 : 



(22) 



2l7lBmT 

In Table [II] the CLFQ results are obtained by first calculating the form factors h(q 2 ),b + {q 2 ) and 
6_(g 2 ) using the covariant light-front approach [47] and k(q 2 ) from the heavy quark symmetry 
relation Eq. (|22p and then converted them into the form-factor set V(q 2 ) and Ao^^iq 2 )- 

Form factors in the CLFQ model are first calculated in the spacelike region and their momentum 
dependence is fitted to a 3-parameter form 

F(0) 



F(q 2 



(23) 



1 — a(q 2 /m 2 B ) + b(q 2 /m 2 3 ) 2 ' 
The parameters a, b and F(0) are first determined in the spacelike region. This parametrization is 
then analytically continued to the timelike region to determine the physical form factors at q 2 > 0. 
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The results are exhibited in Table Hill The momentum dependence of the form factors in the LCSR 
approach can be found in [48], while a slightly different parametrization 

F(0) 



Fi - q ) (1 - q 2 /m 2 B ) [I -a q 2 /m 2 B + b q*/m%] 

is used in the pQCD approach for the calculations of the form-factor q 2 dependence 44]. 

For the calculation in LEET, we have followed [52] to use (±(0) = 0.28 ± 0.04 and C||(0) 
0.22 ± 0.03. For the q 2 dependence, we shall use 

- T ,(0) 



(24) 



Cill(^) 



c 



(1 



j 2 /m 2 B ) 2 



(25) 



For the ISGW2 model, the q dependence of the form factors is governed by Eq. (|A1 



D. Light-cone distribution amplitudes 



The light-cone distribution amplitudes (LCDAs) of the tensor meson are defined as 40j 2 

\ ( e {X) *z a z 13 /-Wv 

f T m 2 T J du e «frW*) P M e ofJJ $J(u) + ' ' ' 



(T(P,\)\qi(yhMx 



o 

,(A> p a\ 

Pu ,„ g v (u) - -z 



(Pz) 2 



Pz 



1 € U Z Z ' ^2 



(Pz) 



(Pzf 



m z T g 3 ( U ) + 0(z z )} , (26) 



(T(P,X)\qi(yh^q2(x)\0) = -if T m 2 T / ' due^+^e^pz* ■P a e*^z s ^- g a 



[u 



(27) 



{T(P,X)\q x (y)a IJa/ q2{x 



M7 (iVF 



m T e aB Z Z 
+ (P W Z„ - P^zJ TTTT^ h t (u) 



+ 



,(A)* a _ (A)* a 



fP* 



r /i3(n) + 0(z 2 



1 (A)* Q £j 

(T(P, AJIft^ga^lO) = -f^ T Jdue^ p y +nP ^ eaP 2 Z pz Z h s (u) , 



(28) 
(29) 



where 53 = 53 + $f - 2^, h t = h t - \{§ T ± + /i 3 ), h 3 = h 3 - and z = y - s. Here 

are twist-2 LCDAs, 3 g v , g a ,ht,h s twist-3 ones, and g 3 ,h 3 twist-4. In the SU(3) limit, due to 



2 The LCDAs of the tensor meson were first studied in [52 

3 Since in the transversity basis, one denotes the corresponding parallel and perpendicular states by Am 
and A±, a better notation for the longitudinal and transverse LCDAs will be &l and $t, respectively, 
rather than $11 and <J>j_. Indeed, the transverse polarization includes both parallel and perpendicular 
polarizations. In the present work we follow the conventional notation for LCDAs. 
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the G-parity of the tensor meson, $T, <? a , h t , h s , g% and /13 are antisymmetric under the 

replacement u — > 1 — u 4C ] . 

Using the QCD equations of motion [3]], l^], the two-parton distribution amplitudes g v ,g a ,h t 
and fi s can be represented in terms of 3>J j_ and three-parton distribution amplitudes. Neglecting 
the three-parton distribution amplitudes containing gluons and terms proportional to light quark 
masses, twist-3 LCDAs g a ,g v , ht and h s are related to twist-2 ones through the Wandzura-Wilczek 
(WW) relations: 

1 



dv 



+ / dv 



(v) 



g™ w {u) = 2u I dv 



+ 2u dv 



hY w {u) = ~(2u-l) 



dv 



dv 



&(v) 



hf w {u) = 3 [u / dv 



*m +u t dv ?m 



(30) 



These WW relations further give us 

1 &T(v) 



1 1 f *ji W 



u ®T (v 



-dv 



v 



K(u) 



V 



dv 



1 



-dv 



dv($i(v) - -h t (v)) = uu 



1 



1 



dv 



1 $^(1;) 



dv 



-dv — u 



v 



1 <&T (v 



-dv 



(31) 



The LCDAs ®7 ± (u,(i) and $ t {u,n) can be expanded as 



*j[ ± («,Ai) = 6u(l-n)^aJ l,±) ' T (/i)Cf 2 (2n-l), 



-,3/2. 



* t («,/i) = 3]Ta^ T ^)P m (2 U -l), 



(32) 



fll _L) T 

where the Gegenbauer moments a) 1 with £ being even vanish in the SU(3) limit, fi is the 
normalization scale and P n (%) are the Legendre polynomials. In the present study the distribution 
amplitudes are normalized to be 

1 1 1 

du(2u- l)*jf(u) = J du(2u- = 1, J^du$ t (u) = 0. (33) 
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Consequently, the first Gegenbauer moments are fixed to be a"'" 1 = a^' T = |. Moreover, we have 



3 J du (2u - 1) g a (u) = J du (2u - 1) g v (u) = 1 , 

o o 
l l 

2 J du(2u - l)h s (u) = Jdu(2u- 1) h t (u) = 1, (34) 



which hold even if the complete leading twist DAs and corrections from the three-parton distribu- 
tion amplitudes containing gluons are included. The asymptotic wave function is therefore 



$JJ(u) = 30u(l - u)(2u - 1), 

and the corresponding expressions for the twist-3 distributions are 

<£» = 5(2u - l) 3 , gf{u) = 10u(l - u){2u - 1) , 
15 



(35) 



and 



hf{u) = — (2u - 1)(1 - 6u + 6u z ) , hf(u) = 15u(l - u)(2u - 1) , 



$f (u) = 5(1 - 6u + 6u 2 ). 



(36) 



(37) 



Note that, contrary to the twist-2 LCDA ± (u), the twist-3 one 3>t(it) is even under the replace- 
ment u — > 1 — u in the SU(3) limit. 

For vector mesons, the general expressions of LCDAs are 



<J>y(x, m) = 6x(l — x) 



l + Y.al V {n)Cll\2x-l) 



n=l 



and 



$ v (x,fj.) = 3 
Likewise, for pseudoscalar mesons, 
$p(x, n) = 6x(l - x) 



2x-l + Y,an' V (v)Pn + i(2x-l) 

n=l 



l + Y,an(»)C?/ 2 (2x-l) 



n=l 



$ p (x,n) = 1. 



(38) 



(39) 



(40) 



E. Helicity projection operators 



In the QCDF calculation, we need to know the helicity projection operators in the momentum 
space. To do this, using the identity 

qKyhsi*) = \{Mq\y)q\x)] + l5 [q 1 (y) l5 q 2 (x)} + Y[qi{yh P q 2 (x)} 

-l P l^\y)l P l^\x)] + \<y pX [q\y)a pX q 2 {x)]} 5a (41) 
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and Eqs. ( M - fl29n . we obtain 



(T(P,X)\f a (y) q l(x)\0) 



due i{uPy+uPx)\ fTm ^ 



(Pzf 



—m T g 3 (u) 



+ 



tfJ,V z 



p 



e u(5 z z 



Pz * (Pz) 2 



f <*,(«) + \^u P al^z p PPz^,^-g a {u) 



1 f± 



+ 9^ K/3 ^ " Kf) rz, 



Pz 

^h(n)+e^z^h s (u) 



' (Pz) 3 
+ 0[{x-yf] 



hAu) 



(42) 



5a 



Since any four momentum can be split into light-cone and transverse components as k^ = 
+ + k^_, we shall assign the momenta 

k% = uEnt + p=n$ + kl , k^ = uPr/ + i< - A:^ , (43) 
AuE AuE 

to the quark and antiquark, respectively, in an energetic light final-state meson with the momentum 

P M and mass m, satisfying the relation P^ = Pn^+m 2 n^/(4P) ~ En^_, where we have defined two 

light-like vectors n± with = (1, 0, 0, —1) and n+ = (1,0, 0, 1) and assumed that the meson moves 

along the direction. To obtain the light-cone projection operator of the meson in the momentum 

space, we take the Fourier transformation of Eq. (|42p and apply the following substitution in the 

calculation 



d 



dki, 



n+ d d 
2E~dvl + dk7 



(44) 



where terms of order k\ have been omitted. The longitudinal projector reads 

2^ 



M,f (A = 0) 



■ St „ 



2P J 



Jt 



r C/xi/ n^n^ h t (u) - iE I dv (&±(v) - h t (v)) a^nt 



dk 



+ ff^fzEa, u n»_eW^5 Xfi 
JT E 

Jt „ 
—i — E 
4 



d 
dk° 



e ^ n+n+ Up; 



St E 



3i f u ( 2 \ 

a^ v n^_n\ht(u) - —E dv \ $±(v) - -^h t {v)\ a^n 



, d h' s (u) 



dk 



(45) 



and the transverse projectors have the form 



f± 

Mj(A = ±1) = -i-^E 



*(A) a (™r\ 
e± ^ n+ { 2E ) 



7^ ji-$l(u) 



+ 



St m T 



^g v (u) - E j U dv (2$J» -g v (vj) fL. 







dk 
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and 



Ml(A = ±2) = -AeI '^Ea^_e^5 K±2 ^- / ~ *,*!(«) + O ( ] L (47) 







E 



dk c , 



(46) 



E 2 



The exactly longitudinal and transverse polarization tensors of the tensor meson, which are 
independent of the coordinate variable z = y — x, have the expressions 



e *(0)^ n H 



fc i "-i. 



2 2£ 2 



3 



c *(A)/w 



m 



AE 2 



1 n_ 



in 



T u 
—n. 



*{\)va + 
t n a U t 
TV 



AE 2 

*(\)ua n 



a u 



(48) 



which in turn imply that 



*(a),„/< t ,- fey 



>A,0 



(1) = VT^ 1 ^' ■ 

The projector on the transverse polarization states in the helicity basis reads 



(49) 



'Ay/2 



El i^liu) 



1 fa mr 
2# E 



^( T l)(l-75) (9v(v)± 



-E^(l- l5 )U dv(2^(v)-g v (v))T 
-E fU(l + 75) ( I" dv(2$T(v) - g v (v)) ± 



+ ^(Tl)(l+75)(ff v (n)T 
3a(«) 



9'a{u) 



<(Tl) 



<9/c 1 „ 



(50) 



After applying the Wandzura-Wilczek relations Eq. (|3ip. the transverse helicity projector ()50p can 
be simplified to 



•4 

4^/2 



£ ^*(Tl) i-<S>l{u) 



fa mr 
ft E 



e*(=Fl)*i(u) 7" (1 T 75) + n£ ^_(1 T 75) 



+ e*(Tl)^(«) 7"(1 ± 75) - uE fl_(l ± 75) 







dk±_ v 



dk\ 



+ 



E 2 



(51) 



to be compared with 

rV r..\ — Jv 



M£(u) 



fv m v 
ft E 



eUTl)*l(u) (>(1 T 75) + uE fi_{l T 7s)^|--) 
+ e*( T l)*^(«) (V (1 ± 75) - n£ ^_(1 ± 75)^|--)] + O ^- 



(52) 
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for the vector meson. The longitudinal projector for the tensor meson can be recast as 



MT (A = 0) 



Jt [2 



4 V 3 



El ft- *f(«) + |^^*t(«) 



i d 



(53) 



to be compared with 



fv 



■■y t t fv m v 



fv E 



2 a nv n^n u + (u — u) - iEuu a^nt 







m 



+ ° # < 54 > 



for the vector meson. 



F. A summary of input parameters 

It is useful to summarize all the input parameters we have used in this work. Some of the input 
quantities are collected in Table HVl 

The Wilson coefficients Cj(/z) at various scales, \i = 4.4 GeV, 2.1 GeV, 1.45 GeV and 1 GeV 
are taken from 54J. For the renormalization scale of the decay amplitude, we choose \x = m&(m;>). 
However, as will be discussed below, the hard spectator and annihilation contributions will be 
evaluated at the hard-collinear scale Hh = y/JIK^ with « 500 MeV. 



III. B -> TP, TV DECAYS 



Within the framework of QCD factorization 24|, the effective Hamiltonian matrix elements are 
written in the form 

{M X M 2 \U^\B) = ^ >^ ) {M l M 2 \T A p ' h +W' h \B) , (55) 

V & p=u,C 

where X p 9 ^ = V p bV* q with q = s,d, and the superscript h denotes the helicity of the final-state me- 
son. For decays involving a pseudoscalar in the final state, h is equivalent to zero. 7~ A p,h describes 
contributions from naive factorization, vertex corrections, penguin contractions and spectator scat- 
tering expressed in terms of the flavor operators a p,h , while l~B p,h contains annihilation topology 
amplitudes characterized by the annihilation operators b^ ,h '. In general, T A p ' h can be expressed in 
terms of ca p i ' h {M l M 2 ) X^ BMl ^ for Mi = T or caf h (M 1 M 2 ) x (BMl ' M2) for M 2 = T, where c 
contains factors arising from flavor structures of final-state mesons, a% are functions of the Wilson 
coefficients (see Eqs. (|B1|) and (|B2|) ). and we have defined the notations 

X (ST,P) = {P \j»\a) {T \J^ = -i2f P A* T (ml)^e*^(0) p BflPBu , (56) 

tub 

X {BP ' T) = -2if T m BPc F* p (m 2 T ), (57) 
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TABLE IV: Input parameters. The values of the scale dependent quantities fy (/i) and a 12 (/u) 
are given for \i = 1 GeV. The values of Gegenbauer moments are taken from [37| and Wolfenstein 
parameters from j^ . 



Light vector mesons [37, 39} 


V 


MMeV) 


/^(MeV) 


a\ 




±,v 
a i 


-L.V 

a 2 


P 
to 

K* 


216 ±3 
187 ±5 
215 ±5 
220 ±5 


165 ±9 
151 ±9 
186 ±9 
185 ± 10 






0.03 ± 0.02 


0.15 ±0.07 
0.15 ±0.07 
0.18 ±0.08 
0.11 ±0.09 






0.04 ± 0.03 


0.14 ±0.06 
0.14 ±0.06 
0.14 ±0.07 
0.10 ±0.08 


Light tensor mesons [401 


T 


/ T (MeV) 


#(MeV) 


||,T ±,T 


/ 2 (1270) 


102 ±6 


117 ±25 


5 
3 


$(1525) 


126 ±4 


65 ± 12 


5 
3 


a 2 (1320) 


107 ±6 


105 ± 21 


5 
3 


iq(1430) 


118 ±5 


77 ± 14 


5 
3 


B mesons 


B 


mB(GeV) 




/s(MeV) 


As(MeV) 


B u 


5.279 


1.638 


210 ± 20 


300 ± 100 


B d 


5.279 


1.525 


210 ± 20 


300 ± 100 


Form factors at q 2 = [37, 41] 


^oT(O) 




A? K *(0) 


^(0) 


V BK *(0) 


0.35 ± 0.04 


0.374 ± 0.033 


0.292 ± 0.028 


0.259 ±0.027 


0.411 ±0.033 


*8T(0) 




Af p (0) 


Af(0) 


C(0) 


0.25 ± 0.03 


0.303 ± 0.029 


0.242 ± 0.023 


0.221 ±0.023 


0.323 ± 0.030 


^oT(O) 




Af^iO) 


A$»(0) 


V Bu (0) 


0.296 ± 0.028 


0.281 ± 0.030 


0.219 ±0.024 


0.198 ±0.023 


0.293 ± 0.029 


Quark masses 


mb(mf,)/GeV 


m c (mb) / GeV 


mP° lc /mP° le 


m s (2.1 GeV)/GeV 


4.2 


0.91 


0.3 


0.095 ± 0.020 


Wolfenstein parameters [381 


A 


A 


P 


V 


7 


0.812 


0.22543 


0.144 


0.342 


(67.2 ±3.9)° 



for the decays B — > TP, and 
*f T,V) = (V\J»\0)(T\J'\B) 



-if v m v 



(4 • e$r)(mjj + m T )Af T (m^) 
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2A BT (m v ) 



2V BT [m\) 



X 



(BV,T) 



m B + rriT tub + tut 

(m% -ml- m T )(m B + m v )A BV (m T ) - ^§-A BV (m T ) 



(58) 



2in\ 



-if T m B m T [(l + ^) A BV (m T 



for h = 0, 



=F 



2p c 



7 a "K) for ft = ±1, 



(59) 



for the decays B — > TV, where x^ 1 ^' and X)^"'^' are expressed in the B rest frame. Note 



ms+mv 

and Xf V > T) 

that in the factorization limit, the factorizable amplitude X^ BMl ' T "> = (T\J fl \0)(M\J'^\B) vanishes 
as the tensor meson cannot be produced through the V — A or tensor current. Nevertheless, 
beyond the factorization approximation, contributions proportional to the decay constant fx of the 
tensor meson defined in Eq. (fT2j) can be produced from vertex, penguin and spectator-scattering 
corrections. 

To evaluate the helicity amplitudes of B — > TV, we work in the rest frame of the B meson and 
assume that the tensor (vector) meson moves along the — z (z) axis. The momenta are thus given 
by 

p T = (E T ,0,0, -p c 



p B = (m B , 0,0,0), 



Pv 



(E v ,0,0, Pc ). 



(60) 



The polarization tensor ej^ of the massive tensor meson with helicity A can be constructed in 
terms of the polarization vectors of a massive vector state 

4"(0) = ( Pc ,0,0,-E T )/m T , e;"(±) = (0,Tl,-2,0)/\/2. (61) 

For the vector meson moving along the z direction, its polarization vectors are 



= (Pc,0,0,E v )/m v , 4 M (±) = (0,=fM,0)/>/2, 



(62) 



where we have followed the Jackson convention, namely, in the B rest frame, one of the vector 
or tensor mesons is moving along the z axis of the coordinate system and the other along the — z 
axis, while the x axes of both daughter particles are parallel 55|. The longitudinal (h = 0) and 
transverse (h = ±1) components of factorization amplitudes X^ ' then have the expressions 

\m\p\ 



X, 



(BT,V) 




ifv 



2m 



TPt 

T 



(m| 



my - m T )(m B + m T )A BT (m v ) 



A BT {m v ) 



X 



(BT,V) 



-if v m B m v 



\f2mj* 



1 + 



m B J 



A BT (m v ) =F 



2 Pc 



m B + mr 



m B + mx 
V BT (m 2 v ) 



Likewise, the factorizable B — > TP amplitude can be simplified to 



X (BT,P) = _ - 2 



2 m B 2 a bt i i \ 
-fp—p 2 c A BT (m 2 P ). 

o TflT 



(63) 



(64) 



The flavor operators a^' h are basically the Wilson coefficients in conjunction with short-distance 
nonfactorizable corrections such as vertex corrections and hard spectator interactions. In general, 
they have the expressions 24j, |57j 



p,h 



(MiM 2 



Ci + 



Cj±l 
N r 



N?(M 2 ) + 



Ci±l CfOs 



N c 4vr 



4tt 



V l h {M 2 ) + — H\ (MiM 2 ) + Pf ,/l (M 2 ), (65) 



N, 



where i = 1, • • • , 10, the upper (lower) signs apply when i is odd (even), q are the Wilson co- 
efficients, C F = (iV c 2 - 1)/{2N C ) with iV c = 3, M 2 is the emitted meson and M\ shares the 
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same spectator quark with the B meson. The quantities (M 2 ) account for vertex corrections, 
Hf{M\M 2 ) for hard spectator interactions with a hard gluon exchange between the emitted meson 
and the spectator quark of the B meson and Pj(M 2 ) for penguin contractions. The expression of 
the quantities iV"/ 1 (M 2 ), which are relevant to the factorization amplitudes, reads 

h, s f° * = 6,8, 
1 else, 



N?(V) 



N?(T) = 0, Ni(P) = 1. 



(66) 



Vertex corrections 



The vertex corrections are given by 



V?(M 2 ) 



C(M 2 ) y dx$j /2 (x) 



C(M 2 ) I dx^ff 2 (x) -121n— + 6-g(l-x) 
Jo 11 L A* 



12 In— -18 + 5 (x) 



1 



I C(M 2 )J 



dx& m2 (x) — 6 + h(x 



V±(M 2 ) 



D(M 2 ) I dx$f 2 {x) 12 In— - 18 + £r(x) 
Jo L A* 



(i = 1-4, 9, 10) 

(t = 5,7), 

(i = 6,8), 
(i = 1-4, 9, 10), 



D(M 2 ) / dx$^ 2 (x) 



12 In— + 6-5r(l-x) (z = 5,7), 
At 

(i = 6,8), 



with 



g{x) 



1 — x 



+ 



2Li 2 (x) -ln 2 x + 



2 In x 
1 — x 



(3 + 2i7r) In x — (x -H- 1 — x) 



and 



/i(x) = 2Li 2 (x) — In 2 x — (1 + 2in) lnx — (x o 1 — x) , 

. . , . lnx 

9t{x) = g{x) + — , 
x 

C(P) = C(F)=Z>(V) = 1, C(T) = y|, D(T) = ^=, D(P) = 0, 



(67) 



(68) 



(69) 



(70) 



3' y/2' 

where x = 1 — x, is a twist-2 light-cone distribution amplitude of the meson M, $ m (for 
the longitudinal component), and &± (for transverse components) are twist-3 ones. Specifically, 
<3? m = 3> t , §> v , & p for M = T,V,P, respectively. The expressions of C(T) and D(T) are obtained 
by comparing Eqs. (f5T) - ([54|) . 



Hard spectator terms 

arise from hard spectator interactions with a hard gluon exchange between the 
emitted meson and the spectator quark of the B meson. Hf(M\M 2 ) have the expressions: 



H KM,M 2) ^^^^^C { M 2 )j\ud V 



uv +r>( C{Mi)uv 



M 2 , 



C(M 1 



(71) 
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for i = 1 - 4, 9, 10, 



C(Mi 



+ r 



$ mi fo)*pfo) 
x C{M 1 )uv 



M 1 



(72) 



for i = 5,7, and J fff(MiM 2 ) = for i = 6,8, where the upper signs are for TV modes and the 
lower ones for TP modes. The transverse hard spectator terms H^{M\M2) read 



M2) _ ^iMkSu,n U2 ^ ^♦J*<.)**<. 



m B X 



(SMi,M 2 ) A S Jo 



F+(M!M 2 ) = 

for « = 1 — 4, 9, 10, and 
iJ-(M!M 2 ) = 

Hf{M 1 M 2 ) = 

for i = 5, 7, and 

H~(M 1 M 2 ) 
H+{M 1 M 2 ) = 0, 



V2if B fM 1 fM 2 inM 1 mM 2 m B f 1 ^ (u - v )$ + 1 2 (u) 



2 X {BM U M 2 ) \ B J Q 



u 2 v 2 



^ifBfM 1 fM 2 mM 1 mM 2 mB_ f l 



U 2 f 



(it-v)^ 1 ^- 2 ^) 



2 X {BM U M 2 ) 



dudv 



V2m B X 



(bm u m 2 ) m 2 M2 X B Jo 



U 2 V 2 



ifBfM 1 fM 2 m M2 m B m Ml m B f 1 ^^(tt)^ 2 ^) 



(73) 
(74) 

(75) 
(76) 

(77) 
(78) 



for z = 6,8. Since we consider only TP and TV modes in the present work, it is obvious that 
M\M 2 = TV or VT for the transverse components. 

Penguin terms 

At order a s , corrections from penguin contractions are present only for i = 4,6. For i = 4 we 
obtain 



G H m 2 (s p ) + 9m 2 
b 



+ C 3 



G h M 2 {s s ) + G h M2 {l) + 2gl 



2cf g G h g , 



4irN c 

+ (c 4 + c 6 )^ [G^ 2 ( Si ) + 5M 2 
where Sj = m 2 /m 2 and the function G^ 2 (s) is given by 

G^(s) = 4 f 1 du<S> M2 > h (u) f 1 dxxx\n[s-uxx-ie], 
Jo Jo 



JM 2 



(79) 



9m 2 



3 /U 3/ Jo 



9'k = t^—[ * m ' h (x)dx, 
6 fl Jo 



(80) 
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with $ M2 '° = C(M 2 )<S>^\ § M2 '± = D(M 2 )<S>± 2 . For i = 6, the result for the penguin contribution 
Cf(x s 



IS 



Pg' h (M 2 



4-irN r 



\ciGm 2 (s p ) + c 3 



+ (c i + c 6 )J2G h M2 (si)\ 



(81) 



In analogy with (|80|) . the function Gm 2 (s) 1S defined as 



G° M Js) 



C(M 2 ) 



du$ ma (u) / dx xx In [s — uxx — ie], 



G^ 2 (s) = 0. 



(82) 



Therefore, the transverse penguin contractions vanish for i = 6, 8: Pg g P = 0. Note that we have 
factored out the r^ 2 term in Eq. f)81 1) so that when the vertex correction V^g is neglected, a[j will 
contribute to the decay amplitude in the product r* f2 a[| 



r M 2 pO 



For i = 8, 10 we find 



a,-, 



9vriV f 



(Cl +iV c C 2 )G^ 2 (s p ;. 



P^(M 2 ) 



Or 



(ci + N c c 2 ) 



G h M (s p ) + 2g 



M 2 



For i = 7, 9, 



7,9 (A*2 



37T 



???. 



Ai 2 



27vr 



777,^ iP* 

5 PC In + 5 PU In + 1 



^ 2 



(83) 
(84) 

(85) 



for M 2 = p , 0J,<f>, and vanish otherwise. Here the first term is an electromagnetic penguin contri- 
bution to the transverse helicity amplitude enhanced by a factor of mBraj/m^, as first pointed 
out in 56]. Note that the quark loop contains an ultraviolet divergence for both transverse and 
longitudinal components which must be subtracted in accordance with the scheme used to define 
the Wilson coefficients. The scale and scheme dependence after subtraction is required to cancel 
the scale and scheme dependence of the electroweak penguin coefficients. Therefore, the scale fi in 
the above equation is the same as the one appearing in the expressions for the penguin corrections, 
e.g. Eq. (f8Q|) . On the other hand, the scale v is referred to the scale of the decay constant fM 2 {v) 
as the operator qj^q has a non- vanishing anomalous dimension in the presence of electromagnetic 
interactions 



57|. The v dependence of Eq. ([85]) is compensated by that of /m 2 ( z/ ) 
The relevant integrals for the dipole operators Si7 are 



<3 



C(M 2 ) 



i $f 2 (u) 
du — 



u 



(86) 



Gf = D{M 2 ) f 1 ^ 
y Jo u 



- ( 



M 2 , 



u$± 2 (u) + - ( u$™ 2 (u) + U$ 



M 2 , 



ai 2 , 



Using Eq. Gf can be further reduced to 



Gl 



g: 



D{M 2 ) 
0. 



du 



^ 2 (u)-^ 2 (u 



M 2 , 



(87) 
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Hence, in Eq. ([H7|) are actually equal to zero. It was first pointed out by Kagan 58J that the 
dipole operators Qs g and Q-j^ do not contribute to the transverse penguin amplitudes at 0(a s ) 
due to angular momentum conservation. 



Annihilation topologies 

The weak annihilation contributions to the decay B — > M1M2 (with M\M.2 = VT or TV) can 
be described in terms of the building blocks b?' h and frfgw 



V * p=u,c V ^ p=«,c i 



(88) 



The building blocks have the expressions 

C F 



bi 

^4,EW 



N 2 



ciA\, 



Cf A i 

m C2Ai > 

c 

Cf 
N 2 
Cp 
N 2 



Cf r 
N 2 
C f 



c 3 A[ + c 5 (4 + A f 3 ) + N c c 6 A 



c 



C 4 ^i + Qj^ 



+ c 7 {A\ + A J 3 ) + N c c 8 Al 
c 10 A\ + c 8 Ai 



(89) 



where for simplicity we have omitted the superscripts p and h in above expressions. The subscripts 
1,2,3 of A 1 ^ denote the annihilation amplitudes induced from (V — A)(V — A), (V — A)(V + A) and 
(S — P)(S + P) operators, respectively, and the superscripts i and / refer to gluon emission from 
the initial and final-state quarks, respectively. Following 57|] we choose the convention that M2 
contains an antiquark from the weak vertex with longitudinal fraction v, while M\ contains a quark 
from the weak vertex with momentum fraction u. The explicit expressions of weak annihilation 
amplitudes are: 



4°(MiM 2 ) 



-vra s ^ dudvl ^ h (u)^ h (v) 



1 



+ 



1 



T ^r^^ mi (u)^ m2 (v) 



u(l — uv) uv 2 
2 ) 



uv 



4l'-(MiM 2 ) 
Ai' + tMiMa) 



-vra. s 



in 



B 



dudvl ^{u)^ 2 {v) 



M 2 . 



U + V 



+ 



u 2 v 2 (1 — uv)- 1 



-7ra, 



■m Ml mM 2 



in 



B 



dudvl 



M 2l 



uv 3 (1 — uv) 2 v 2 {\ — uv) 



A i i°{M 1 M 2 ) 



-ira s f dudvl^iu)^ 12 ^) 
3 Jo 



1 1 

+ 



v (1 — uv) ' u 2 v _ 
2 1 



(93) 



A\ (MiM 2 ) = -vnx 



\/2m Ml mM 2 



B 



J^dudv\^{u)^ 2 {v) 
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U + V 1 



v?v 2 (1 — uv) 2 



4' + (M!M 2 ) = -7m s ^ mM 2 imM2 / ,1 d U (fo(*^(«)^ a («) 



u 3 u (1 — mi;) 2 u 2 (l — uv) 



(94) 



(95) 



+ 



C(M 



W(l — TO J 



(96) 



4' - (MiM 2 ) 



./<> 



m w , x m>(l 



m M2 



uv 



uv(l — uv) 



(97) 



A^°{M 1 M 2 ) 



7ra, / dudv 
h 



\C(M 1 
C(Mi 



2J_ Mi 

x 



C(M 2 



r M 2 ^Mi 

X 



$1 («)* 




U 2 t> 



(98) 



i^' _ (MiM 2 ) 



and ^{'^ = A 



lf,h 



A 



ira. 



V2 



du dv 



m M , 



2 ~Ml fT> A/ l 



-r 



x 



+ 



rf 2 ^ Jl (u) $f» 



m M2 



4 



(99) 



= 0. Here in the helicity amplitudes with h = 0, the upper signs 
(T,V), (V,T), and (V,P) and the lower ones to (M 1 ,M 2 ) = (P,V). 



correspond to (Mi,M 2 ) 

When (Mi,M 2 ) = (V,P), one has to add an overall minus sign to A % f . For (Mi,M 2 ) = (P,V), 
one has to change the sign of the second term of A^ . Note that in this paper, we adopt the 



(*./).o 



for the TP modes. 



notations A 

j j 

Since the annihilation contributions A{ 2 are suppressed by a factor of mim 2 /m^ relative to 
other terms, in the numerical analysis we will consider only the annihilation contributions due to 



A f,o A f,- 



^1,2,3 anc ^ -^3 



Finally, two remarks are in order: (i) Although the parameters ai{i ^ 6,8) and a%fir x are 
formally renormalization scale and 75 scheme independent, in practice there exists some residual 
scale dependence in aj(/i) to finite order. To be specific, we shall evaluate the vertex corrections 
to the decay amplitude at the scale \x = m^. In contrast, as stressed in [24|, the hard spectator 
and annihilation contributions should be evaluated at the hard-collinear scale \ih = yji^h with 
Ah ~ 500 MeV. (ii) Power corrections in QCDF always involve troublesome endpoint divergences. 
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For example, the annihilation amplitude has endpoint divergences even at twist- 2 level and the hard 
spectator scattering diagram at twist-3 order is power suppressed and possesses soft and collinear 
divergences arising from the soft spectator quark. Since the treatment of endpoint divergences is 
model dependent, subleading power corrections generally can be studied only in a phenomenological 
way. We shall follow j^] to model the endpoint divergence X = Jq 1 dx/x in the annihilation and 
hard spectator scattering diagrams as 

X A = In (^j (1 + PA e^), X H = In (1 + p H e^), (100) 

with the unknown real parameters pa,h an d <\>a,h- F° r simplicity, we shall assume that X\ and 
Xfl are helicity independent; that is, X^ = Xj[ = X A and Xjj = X^ = Xjj. 



IV. NUMERICAL RESULTS 



Let the general amplitude of B — > TP be 



A^ TP = ^{aX^ + -aX^). (101) 



Its decay rate is given by 

i\ 

'<ij 

It follows from Eqs. ([57]) and that 



r _ _ Gjp_ fppl 

1 B-s-TP 



6tt mji 



aA Q (m P ) + y-a 



BP 1^2 



(103) 



2 f P p c A^ T (ml) 

where p c is the center-of-mass momentum of the final-state particle T or P. Note that the coefficient 
a vanishes in naive factorization. 

The decay amplitude of B — > TV can be decomposed into three components, one for each 
helicity of the final state: Ao,A+,A-. The transverse amplitudes defined in the transversity basis 
are related to the helicity ones via 

The decay rate can be expressed in terms of these amplitudes as 

B B 
Writing the general helicity amplitudes as 

A ± (B -+ TV) = ^ (b±xi B Ty) + V ' T) ) , (107) 

where Xq ± ' and X ± ' are given in Eqs. (|63|) and ([59]) . respectively, and it is understood that 
the relevant CKM factors should be put back by the end of calculations, the decay rate has the 
following explicit expression 



Ao(B -> TV) = % ( b°X^ + Wxf V ' T) ) , (106) 



T B->TV = jf + fol + 7Pc + Ap c a ), (108) 

487T nij, 
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with 



mr, 



a 

P = 6 
7 = 6 
A 



2 2 

myrrij, 
(tub + bit) 



2 l w A 2 I ' 



&+ |2 (v BT ) 2 + | 6 -|2 (v BT ) 2 



4(m| 



4 



|6°| 2 Af T Af T , 



TUB ' 



(109) 



(me + mr) 



2m| 



(Af T ^ 2 



+ |&-| 2 (Af 7 ;) 2 )m 2 / 774 + |& u | 2 (m 



4 



4 



m|) 2 (A 1 



BT\2 



where we have adopted the shorthand notations, 



Af T 



A BT 



BT 
L l,± 



^f T (m 2 y 



+ 



b° [3 f T mr m T m B + my ^ BV 



77> 



A; 



b° y 2 fv my p c mfl + my 
U V 2 /y my p c mg + my 



A? T (ml] 



6 ± /y my p c m B + m T 
b ± 



(110) 
(111) 
(112) 
(113) 



b ± " fv my p c tub + my 

Note that Eqs. (I103P and (I108P are in agreement with [161 ] for the special case that a = b° = b ± = 1 
and a = W = b^ = 0. As stressed in [16j], the ^ dependence in Eq. (|103p indicates that only the 
L = 2 wave is allowed for the TP system, while in the TV modes the L = 1, 2 and 3 waves are 
simultaneously allowed, as expected. 



A. B -> PT decays 



As noticed in [59] , since the penguin annihilation effects are different for B — > VP and B — > PV 
decays, the penguin annihilation parameters X\ p and X PV are not necessarily the same. Indeed, 
a fit to the B -> VP,PV decays yields ps 1.07, < p « -70° and p£ v » 0.87, ^ y « -30° 



59J. Likewise, for B u d —¥ TP decays we find that the data of B u ^ — > TP can be described by the 



penguin annihilation parameters 

P Y = 0.83, f/ = -70°, p7 = 0.75, 4> p a t = -30° . (114) 

For B — )• T transition form factors, the LEET or pQCD predictions are favored by the exper- 
imental data of B -»• / 2 (1270)if and B -)• / 2 (1270)vr, while the CLFQ or ISGW2 model results 
are preferred by the measurements of B — > K|(1430)?/'\ (1430)w, -ftT^ (1430)</>. For example, the 
branching fractions (in units of 10 -6 ; only the central values are quoted here) for the / 2 vr~, / 2 ^° 
and feK~ modes are found to be 8.1 (2.7), 5.0 (3.4) and 6.4 (3.8), respectively, using the CLFQ 
(LEET) model for B — )■ / 2 transition form factors. The corresponding experimental values are 
1.57^49, 2.71;};!, and i- 06 ^!!- Therefore, it is evident that the data favor LEET over the CLFQ 
model for decays involving B — > / 2 transitions. Likewise, the branching fractions for the modes 
K^cj), K^rj and K 2 rf are calculated to be 7.4 (4.7), 6.8 (4.7) and 12.4 (8.4), respectively, using 
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TABLE V: CP-averaged branching fractions (in units of 10~ 6 ) and direct CP asymmetries (%) 
for B — > PT decays with AS" = 1. The parameters pa and (ft a are taken from Eq. (|114p . The 
theoretical errors correspond to the uncertainties due to the variation of Gegenbauer moments, 
decay constants, quark masses, form factors, the Xb parameter for the B meson wave function and 
the power-correction parameters pa,h, 4>A,H- Then they are added in quadrature. The experimental 
data are taken from j6o| and the model predictions of (20I are for = 0.3 . 

Decay — Acp 

QCDF KLO [20] MQ [21] Expt. 

B~ -+ if 2 *(1430)°7r- 3.1±|;f 5.6±f;| L6+fj 

B~ -> K^(U30)-tt° 2.2+f;^ 0.090 0.15 0.2±J£f 

B° -> i^ 2 *(1430)-vr+ 3.3±1;| < 6.3 1.7±f J 

B° -> if 2 *(1430)°vr 1.2±f;| 0.084 0.13 < 4.0 7.1±H;f 

B~ -> a 2 (1320)°K- 4.9±|;| 0.311 0.39 < 45 -27.lt|?;f 

5~ -> a 2 (1320)-K° 8.4^ 0.011 0.015 -0.6i°;8 

5° -»■ a 2 (1320)+if- V-^H 0.584 0.73 — 21.5±||;g 

B° a 2 (1320)°K° 4.2±|;| 0.005 0.014 6.7±|;| 

B~ -»■ / 2 (1270)i^- 3.8±^ 0.344 l-06l{j:ig -39.5±|§;| 

5° -> / 2 (1270)i?° 3.4±l;f 0.005 2.7±}; 2 1 — 7.3±f;| 

B~ -> /^(1525)A- 4.0j£| 0.004 < 7.7 -O^o 

B° -> /^(1525)F° 3.8l^| 7 x 10~ 5 0.8l^ 

B~ -»■ if|(1430)-r? SS- 1 ^ 0.031 1.19 9.1 ±3.0 l-Sl^e 

B~ -»■ if|(1430)-»/ 12.1±?g;J 1.405 2.70 28.0±|;g -1.7±|;| 

B° ->■ X|(1430)°77 6.6± 1 | f 0.029 1.09 9.6 ±2.1 S-^ti 

B° -> X|(1430)V 12-4t^j 1.304 2.46 13.7±|;? — 2.2±|;g 



the CLFQ (LEET) model for B — >• i^l transition form factors. The corresponding experimental 
values are 8.4 ± 2.1, 9.1 ± 3.0, and 13.7^!'^. It is clear that the CLFQ model works better for 
decays involving B — > transitions. In this work we shall use the B — > if|(1430) form factors 
obtained in the CLFQ model and B — > a 2 (1270) and B — > / 2 ones from LEET (see Table [TT]) . 

Branching fractions and CP asymmetries for B — > TP decays are shown in Tables [V] and IVIl 
The theoretical errors correspond to the uncertainties due to the variation of (i) the Gegenbauer 
moments, the decay constants, (ii) the heavy-to-light form factors and the strange quark mass, 
(hi) the wave function of the B meson characterized by the parameter As, and (iv) the power 
corrections due to weak annihilation and hard spectator interactions described by the parameters 
PA,H, 4>A,H- We allow the variation of pa and (ft a to be ±0.4 and ±50°, respectively, and put pu 
and (ftjj in the respective ranges < pu < 1 and < (ftjj < 2-7T. To obtain the errors shown in these 
tables, we first scan randomly the points in the allowed ranges of the above-mentioned parameters 
and then add errors in quadrature. Power corrections beyond the heavy quark limit generally give 
the major theoretical uncertainties. 

For B —7- K\rf- ) decays, there exist three different penguin contributions as depicted in Fig. [TJ 
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TABLE VI: Same as Table |V] except for B -+ PT decays with AS = 0. 
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u - ±z -0.12 
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5°^ 


K$(143Q)+K- 
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(i) b — > sqq — >■ sr) q , (ii) b — >■ sss — > srj s , and (iii) b — > sqq — > qK%, corresponding to Figs. 1(a), 
1(b) and 1(c), respectively. The dominant contributions come from Figs. 1(b) and 1(c). Since 
the relative sign of the r] s state with respect to the r\ q is negative for the rj and positive for the 
rj' (see Eq. (|D1|) ). it is evident that the interference between Figs. 1(b) and 1(c) is destructive 
for K^rj and constructive for K^r/' . This explains why K^rj' has a rate larger than K^rj. It was 
known that the predicted rates in naive factorization are too small by one order of magnitude, of 



order 1.0 x 10~ 6 for B(B° -+ Kfi]) and 2.5 x 10~ 6 for B(~B° -+ Kfrf) [2lJ, 4 One reason is 
that the factorizable contribution to Fig. 1(c) vanishes in the naive factorization approach. The 
rates of K^vf^ are greatly enhanced in QCDF owing to the large power corrections from penguin 
annihilation and the sizable nonfactorizable contributions to Fig. 1(c). 



4 The rate of 75° — > K^r\ was predicted to be very suppressed in [20| (see Table |VJ due to the use of a 
wrong matrix element for (f/ ^|s75s|0) [sij - 
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TABLE VII: CP-averaged branching fractions (in units of 10 -6 ), direct CP asymmetries (%) and 
the longitudinal polarization fractions /l for B — > VT decays with AS = 1. The parameters pa 
and (ft a are taken from Eq. (|124j) . 



Decay 












JL 


^CP 




KT,0 [901 














OCDF 


MO [911 


Fvnt 


OCDF 


T7 Y nt 


JD — ) 


A 2 (143U) p 


18.0_ 17 _2 








n fi^+o io 

u - D,i -0.09 




1 o+ - 8 

J - u -1.0 


D — ) 


A 2 (14,30) p 


1U.4_ 9 2 




n 74 




u - DD -0.07 




n 1+11.1 
z - i - 9.9 


Tf _v 

D — 7 


A ic\\— „+ 
A 2 (,14*51) J P 


lq c+52.0 








fi4+ - 07 
u - D4 -0.03 




-1 5+ 2 - 6 


if -A 

D — r 


A 2 ^14oUJ p 


q 5+33.4 
y -°- 9.5 


D 9^ 


n R8 




64+ - 15 

u - D ^-0.37 




-4 0+ 14 - 1 

4 - u -10.8 


D— v 

D — 1 


i\ 2 ^14oUJ 


7 r: + 19.7 
7.0 


n 1 1 9 


u.uu 




64+ - 08 
u - D4 -0.07 


n k« _i_ n 1 1 

U.OU I U.ll 


9 n+12-2 

z - u -10.5 


D — 7 


A 2 ^14oUJ OJ 


8 1+ 2 F 

7.6 






10.1 db 2.3 


66 +0 U 

u -°°-0.15 


n 4.5 4- n -1 

L/.^rtJ —I— Uil^ 


4 4+10-9 


-D — 7 


1\ 2 ^l^toU J (p 


7 ,1+25.8 
5.2 


9 180 


Q 94 


8 4 4-91 

O . 4: _1_ _■ . ± 


85 +ai 2 


n on 1 n 1 n 

U.OU _l_ U.1U 


1 +L2 
u -- L -0.5 


B° —7- 


K|(1430)V 


' •' - 5.5 


2.024 


8.51 


7.5 ± 1.0 


86+°,-^ 
u.ou_ 77 


901+^ 

u - aui -0.069 


09 + S'o? 

U.U»_q 21 


B~ -» 


a 2 (1320)°A"*- 




1.852 


2.80 




n 70+0.22 
u - '°-0.33 




-1 5 0+ 56 - 


B~ -> 


a 2 (1320)-K*° 


c 1+23.8 
D1 - 5.4 


4.495 


8.62 




n 79+0.20 

u - ' y -0.64 




-0 1 +L3 
u - x -0.3 




02(1320)+^*- 


6-l +2 5 4 3 3 


3.477 


7.25 








-13.3l 3 ^o 




a 2 (1320)°K*° 


3.4lf| 


2.109 


4.03 




0.82^ 




1 2+ 7: ° 




f 2 (l270)K*- 


O Q + 17.3 

°- 6 - 6.7 


2.032 






o.93t^ 3 7 




-S.ll 1 , 3 ! 




/ 2 (1270)iT° 


Q 1+18.8 
y - X - 7.3 


2.314 






n n 4 +0.06 
U - y4 -0.69 




-0.08+f 3 


B- -> 


fi(1525)K*~ 


12.6±?f;? 


0.025 






U - DO -0.38 




o.eill 




/^(1525)F* 


ia c:+25.4 
10.0_ 11 g 


0.029 






fifi +0 - 27 
u - DD -0.38 




0.2l°; 3 



b Q 1 b rvf^Ho s 6 rJ^vVu s 

_^ ^2 c _^ ^2 < Vq 

(a) (6) (c) 

FIG. 1: Three different penguin contributions to B — > K^rf \ Fig. 1(a) is induced by the penguin 
operators 03,5,7,9- 



From Tables |V] and IVII we see that the predicted branching fractions for penguin-dominated 
B — 7 TP decays in QCDF are larger than those of 20J and 21[ by one to two orders of magnitude 
through the aforementioned two mechanisms for enhancement, while the predicted rates in QCDF 
are consistent with 20] for the leading tree-dominated modes such as a 2 7r~, a^vr - , /2vr _ . Note that 



K2 7r + vanish in naive factorization, while 



the branching fractions of B~ — > iCJ v ir~ and B 
experimentally it is (5.6l 2 ' 2 ) x 10 -6 for the former. The QCDF calculation indicates that the 
nonfactorizable contributions arising from vertex, penguin and spectator corrections are sizable to 
account for the data. 
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TABLE VIII: Predicted branching fractions (in units of 10 6 ), direct CP asymmetries (%) and the 
longitudinal polarization fractions fi for B — > VT decays with AS" = 0. 
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Just as B —> oott with ao = ao(980) or ao(1450) and S — )• &i(1235)7r decays, we see from Table 
IVII that for B — > a2(1320)7r decays, the and ci^ft modes are highly suppressed relative to 

ci2 7T~ and a^-K~ , respectively. Since (see Appendix B) 

A Bo^n+ K Yl VpbV* d [5 pu ai + al]X {B * A2 \ 

p=u,c 

Ab^4k- « E V pb V P*d [S P u ax + a\\ X^ a ^\ (115) 

p=u,c 

it is tempting to argue that T(B° — > a^ 71 " ) ^ r(B° — >• a^vr + ) is a natural consequence of 
naive factorization as the tensor meson cannot be created from the V — A current. However, the 
suppression of vr + relative to a^ 71 " - ^ n QCDF stems from a different reasoning. The amplitude 
^(Bn,a2) ^ Qes no £ ygjjjgij [ n QCDF owing to the nonfactorizable corrections. Indeed, x^ B7T ' a ' 2 ^ = 
0.80 and X^,*) 

= 0.69 are numerically comparable. Therefore, one may wonder how to see the 
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aforementioned suppression ? The key is the quantity iVj(M2) appearing in the expression for the 
effective parameter Oj [see Eq. (|65p ]. This quantity vanishes for the tensor meson [cf Eq. (|66p ]. As 
a result, the parameter 01(7702) is not of order unity as it receives contributions only from vertex 
corrections and hard spectator interactions, both suppressed by factors of a s /(47r). Numerically, 
we have 01(7702) = —0.035 + i0.014. By contrast, ai(a27r) is of order unity. This explains why 
B° —7- a^Tr - has a rate greater than a2~7r + and why B~ — > a 2 ir° is suppressed relative to o^tt - . 5 
The same pattern also occurs in B — )■ 02/5 decays, see Table IVlTTl 

The branch ing fractions of B — > a^^'-' of order 10~ 7 in QCDF are in sharp contrast to the 
predictions of 2l|], ranging from 25 x 10~ 6 to 70 x 10~ 6 (see Table |VT|) . It seems to us that it is 
extremely unlikely that the rate of a 2 r]( ) can be greater than a2~7r by four orders of magnitude 
as claimed in (2l| . It appears that the former should be slightly smaller than the latter in rates. 
This can be tested in the future. It is also interesting to notice that while B — > K 2 K decays 
are very suppressed in naive factorization, their branching fractions are a few xl0~ 7 in QCDF. 
Finally, it is worth remarking that B° — > K 2 + K~ and B° — > K 2 ~ K + can only proceed through 
weak annihilation. 



B. B ^TV decays 

Branching fractions, direct CP asymmetries and the longitudinal polarization fractions for B — > 
TV decays are shown in Tables IVIII and IVIII1 Thus far only four of the B — > TV decays have 



been measured: B — > K 2 (<f>, uj) and B 



J, co). They can be used to fix the penguin 



annihilation parameters. From Eqs. ()B38|) and (|B40p we have 
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(116) 



Since X [ Bul ' K2) / xf K * 2 ^ ] = 0.56, it is expected that B(B~ K*-oj)/B(B~ -> K* 2 



\X 



(Bu>,K 2 ) /Y (BK 2 , 
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0.15 , provided that the penguin-annihilation parameters are the same 
i.e. 9/?^ = P^a % an< ^ likewise for ^4. However, it is the other way around 
experimentally: the rate of K 2 u) is larger than that of K 2 §. Since in the B — > VV sector, 



for K 2 uj and K 2 



B(B~ —7- K*~oj)/B(B~ —7- K*~(f>) ~ 0.3 60J, it is thus puzzling as why K 2 uj behaves so differently 
from K*u in terms of branching fractions. It is clear from Eq. (|116p that the B — > K 2 <p decay 
receives penguin annihilation via and 4>a ■> while B — > K 2 oj is governed by p^ T and (j)\ T ■ 
Therefore, we should have p\ T p^ in order to account for their rates (see Eq. (|124p below). 

The branching fractions of the tree-dominated modes 02^, /2</>, fi§ ar e very small, of order 
10 -9 (see Table [VTTTT) . as they proceed only through QCD and electroweak penguins. 

For charmless B — > TV decays, it is naively expected that the helicity amplitudes Ah (helicities 



-ft" 



The same argument also explains the suppression of B — > b 1 7r + relative to b^ir in QCDF [621 ] . 
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h = 0, — , + ) for both tree- and penguin-dominated B — > TV decays respect the hierarchy pattern 

Aqcd \ _ / Aqcd 



A : A- : A+ = 1 : : . (117) 

\ m h J \ m b ) 

Hence, they are dominated by the longitudinal polarization states and satisfy the scaling law, 
namely [581 ] . 

f T ^l-h = o(^f], L = 1 + (^), (118) 
V m B J f\\ V m B J 

with /|| and /t being the longitudinal, perpendicular, parallel and transverse polarization 

fractions, respectively, defined as 

fa ^T = iai 2 + i^i 2 + i.4±i 2 ' (119) 

with a = L, ||, _L. 

The so-called polarization puzzle in B — > VV decays is the enigma of why the transverse 
polarization fraction fa in the penguin-dominated channels such as B — > (j)K* , pK* is comparable 
to /l, namely, /t//l ~ 1- This poses an interesting challenge for any theoretical interpretation. 
For B — > TV decays, the experimental measurement indicates that fa/ ft 1 for B — > </>A|(1430), 
whereas Jt/Il ~ 1 f° r B — > ojK^ (1430), even though both are penguin-dominated. 

Consider the ratio of negative- and longitudinal-helicity amplitudes 

A_ 

A 



b-^k;~ 

A_ 

Aq 

The longitudinal polarization fraction can be approximated as 




(120) 



\a c f +/3 3 -| 2 |X^| 2 

Eh=o,-l«f + /3 3 TI^| 

|«r +/3 3 -| 2 |X7 ' 2 



/l(-RT 2 » ~ 1 — ^ ^±1 . (121) 



We have 



: \X uK *\ : |A^*| = 1 : 0.51 : 0.04, 
\X°K*J ■■ \X«*J : |A+^| = 1 : 0.38 : 0.06 . (122) 

In the absence of penguin annihilation, we find fa^uK^) ~ f^K^) ~ 0.72. As we have stressed 



m 



63J, in the presence of NLO nonfactorizable corrections e.g. vertex, penguin and hard spectator 



scattering contributions, the parameters a- 1 are helicity dependent. Although the factorizable 
helicity amplitudes A , A~ and A + or A°, A , A + respect the scaling law (|117p with Aqcd/"i& 
replaced by 2mv ) T/ m B for the tensor and vector meson productions, one needs to consider the 
effects of helicity-dependent Wilson coefficients: A-/Ao = f(a^)X~ /[f(a^)X ]. The constructive 
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(destructive) interference in the negative-helicity (longitudinal-helicity) amplitude of the penguin- 
dominated B — > TV decay will render f(aj~) 3> /(a- 1 ) so that A~ is comparable to Aq and 
the transverse polarization is enhanced. Indeed, we find /t^ujK^) — fri^K^) ~ 0.28 . Therefore, 
when NLO effects are turned on, their corrections on a~ will render the negative helicity amplitude 
A-{B — > K%(p) comparable to the longitudinal one Ao(B — >• K^fj)) so that even at the short-distance 
level, fi for B° — > K^cj) is reduced to the level of 70% and likewise for B° — > K^lo. 

As noticed in passing, penguin annihilation is needed in order to account for the observed rates. 
This is because, in the absence of power corrections, QCDF predicts too small rates for penguin- 
dominated B — > TV and VV, V A decays. For example, the calculated B — > KZ,4> rate is too small 
by a factor of 2.5 and B — > K%uj by two orders of magnitude. We shall rely on power corrections 
from penguin annihilation to enhance the rates. A nice feature of the (S — P)(S + P) penguin 



annihilation is that it contributes to Aq and A- with the same order of magnitude 64 ] 



A P A : -4 PA :A P + A =( ^ln H*V : ( ^ In V : ( **™) * . (123) 
The logarithmic divergences are associated with the limit in which both the s and s quarks orig- 



inating from the gluon are soft 64J. As for the power counting, the annihilation topology is of 
order l/m b and each remaining factor of l/m& is associated with a quark helicity flip. The fact 
that A^ and Aq A have the same power counting explains why penguin annihilation is helpful to 
resolve the polarization puzzle. The relative size of A^_ A and Aq A depends mainly on the phase 
4>A- It turns out that the longitudinal polarization fraction for B — > K%(j) is quite sensitive to the 
phase 4>a V , while fii^K^) is not so sensitive to <^ T . For example, = 0.88, 0.72, 0.48, 

respectively, for 4> T / = -30°, -45°, -60° and /i(wf 2 *) = 0.68, 0.66, 0.64, respectively, for 
= —30°, —45°, —60°. Hence, we can use the experimental measurements of /x to fix the 
phases <p\ T and 0^ and branching fractions to pin down the parameters p^iVT) and pa(TV): 

p T / = 0.65, <$ V = -33°, p^ T = 1.20, <^ T = -60°. (124) 

It should be stressed that although the experimental observation of the longitudinal polarization 
in B — > K^cf) and B — > K^ui decays can be accommodated in the QCDF approach, no dynamical 
explanation is offered for the smallness of fxiK^cft) and the sizable /t(wX|). 

For penguin-dominated B — > TV decays, we find /i(if|p) ~ /i(if|^) ~ 0.65, whereas 
/zCfe-K"*) ~ 0.93 (cf. Table IVIip . It will be very interesting to measure Jl for these modes 
to test the approach of QCDF. Theoretically, transverse polarization is expected to be small in 
tree-dominated ~B — > TV decays except for the a<7p°, p + , K£°K*~ and K%°K* modes. 



V. CONCLUSIONS 

We have studied in this work the charmless hadronic B decays with a tensor meson in the final 
state within the framework of QCD factorization. Due to the G-parity of the tensor meson, both 
the chiral-even and chiral-odd two-parton LCDAs of the tensor meson are antisymmetric under 
the interchange of momentum fractions of the quark and anti-quark in the SU(3) limit. The main 
results of this work are as follows: 
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We have worked out the longitudinal and transverse helicity projection operators for the 
tensor meson. They are very similar to the projectors for the vector meson. Consequently, 
the nonfactorizable contributions such as vertex, penguin and hard spectator corrections to 
B — > T(P, V) decays can be directly obtained from B — > VP, VV ones by making some 
suitable replacement. 

The factorizable amplitude with a tensor meson emitted vanishes under the factorization 
hypothesis owing to the fact that a tensor meson cannot be created from the local V — A and 
tensor currents. As a result, B~ — > K^ti~ and B° — > K^tt^ vanish in naive factorization. 
The experimental observation of the former implies the importance of nonfactorizable effects. 

Five different models for B — > T transition form factors were considered. While the pre- 
dictions of B — > /2(1270) form factors based on large energy effective theory or pQCD 
are favored by experiment, the covariant light-front quark model or the ISGW2 model for 
B — > KK1430) ones is preferred by the data. 

For penguin-dominated B — > TP and TV decays, the predicted rates in naive factorization 
are normally too small by one to two orders of magnitude. In QCDF, they are enhanced by 
the power corrections from penguin annihilation and nonfactorizable contributions. 

There exist three distinct types of penguin contributions to B — > K^i]^'- (i) b — > sqq — > srj q , 
(ii) b — > sss — > si] s , and (hi) b — > sqq — > qK^ with r, q = {uu + dd)/\/2 and t] s = ss. The 
dominant effects arise from the last two penguin contributions. The interference, constructive 
for K^rj and destructive for K^rj between type-(ii) and type-(iii) diagrams, explains why 
F(B -> K$rf) » T(B -> Kir,). 

We use the measured rates of K%oj and K^cj) modes to extract the penguin annihilation 
parameters p T ^ and p^ T and the observed longitudinal polarization fractions f^K^uj) and 
/i(i^ 2 * ( / ) ) to fix the phases <j)^ T and <\?a ■ The unexpectedly large rate of B — > K%ou relative 
to B — >■ K^cj) implies that p^ T 3> p^ T . However, it may be hard to offer more intuitive 
understanding for the large disparity between p T ^ and p^ T in magnitude. 

The experimental observation that /t/Il "C 1 for B — > (pK^i 1430), whereas fr/fi ~ 1 for 
B — > a;-R'2*(1430), can be accommodated in QCDF, but cannot be dynamically explained at 
first place. For penguin-dominated B — > TV decays, we find fiiK^p) ~ /i(i^|a;) ~ 0.65 and 
fh(K*f2) ~ 0.93. It will be of great interest to measure fi for these modes to test QCDF. 
Theoretically, transverse polarization is expected to be small in tree-dominated B — > TV 
decays except for the a^p , a^p + , K^K*^ and K^K* modes. 

For tree-dominated decays, their rates are usually very small except for the 
a^in - , p~), (tt~ , p~) and f2(^~,P~) modes with branching fractions of order 10~ 6 or 
even bigger. 
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Appendix A: Form factors in the ISGW2 quark model 



Consider the transition B — > T in the ISGW2 quark model 4a ]. where the tensor meson T has 
the quark content q\q2 with q~2 being the spectator quark. We begin with the definition 46] 
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/m T \ 1/2 
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18 v 



16 



j bt uiBrriT \33 — 2rij 

rh is the sum of the meson's constituent quarks' masses, m is the hyperfine-averaged mass (for 
example, tub = jThb* + \nriB), t m = (thb — wt) 2 is the maximum momentum transfer, and 
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with ni\ and ni2 being the masses of the quarks q\ and q~2, respectively. In Eq. (|Aip . the values of 
the parameters (5b and (3t are available in [46J] and fi BT = \{(5 B + Pt)- 

The form factors defined by Eq. (|17p have the following expressions in the ISGW2 model: 
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and 



- _ 1 = *m (A6) 

In the original version of the ISGW model [45], the function F n has a different expression in its 
{t m — t) dependence: 




exp 
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where k = 0.7 is the relativistic correction factor. The form factors are then given by 
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Note that the expressions in Eq. (|A4p in the ISGW2 model allow one to determine the form factor 
6_, which vanishes in the ISGW model. 



Appendix B: Decay amplitudes 



The coefficients of the flavor operators a?'^ can be expressed in terms of a\ in the following: 

ai(MiM 2 ) = a x (M 1 M 2 ), 
a 2 (M x M 2 ) = a 2 (M 1 M 2 ), 
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al (AfiAf 2 ) - a p 5 (MiM 2 ) for M X M 2 = TP , 
of (MiM 2 ) + af {MM) for MiM 2 = FT , 
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for B — > TV decays with (M1M2) = (TV) or (VT). It should be noted that the order of the 

X 



arguments of oS(M\M2) and a?(MiM2) is relevant. The chiral factor r^f is given by 



rZ' K ^)= , w 2m '' X WA , r$ l (/i) = -~ ^ (B3) 

x m b {n) (m qi + m Q2 ){fi) f { q %m b {^)m q ^) 



for the pseudoscalar mesons, 

r fiA = 
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= ^- f -M- (B4) 



for the vector meson, and 

„T(,a - 2l7lT /tM 



r T^ ) = ^_JT^± (B5) 



for the tensor meson. See Appendix [D] for further discussions on the parameters hq iS and fq} s - 

In the following decay amplitudes, the order of the arguments of c^ ,h \M\M2) and /3f (M1M2) 
is consistent with the order of the arguments of xiR Ml,M ^ or X^ Ml ' M2 \ where 

Pf(TP) = ~ tf JL fTfp i$, (3f(PT)= ~ lf ^ fTfp bl for TP modes, (B6) 
Pf h (TV) = lI^Mv^h pP,h? VT) = ifBfrfv ph for TV modeg _ 

Ht K ' V (BT,V) 1 ' H% v ; WBV,T) * ' 

The decay amplitudes for P — > TP, TV are summarized as follows: 
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the amplitudes for B — > $21)' can be obtained from B — > jif] with the replacement (fi,rj) ~~ ^ (fiiV')i 
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"yl E ^ | a 4 2 a4 ' EW + ^3 + /^4 2^3,EW 2 P4 ' EW 
-icfBfK'fK 

G F 



(B16) 



(B17) 

(B18) 
(B19) 



(B20) 



(B21) 



x {bk* 2 ,k) 



4,EW 



(B22) 



rr E 4*°' 
v2 p=U}C 

-icf B fKzfK 



«4 - o<EW + 03 + ^4 - o^3 



2^3,EW 2^ 4 ' EW 



(BK,K* 2 ) 



4,EW 



(B23) 



6. Decay amplitudes with AS = 1 



B~^f 2 K- 



Gf_ 

V2 



E a?< 



p=u,c 



Spu «2 + 2af + -or 



p 

3,EW 



■^{BKJD 
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+ 



<5 pu (ai + h) + <4 + < EW + /f + /3| EW ] X^-*), 



V2A 



B°^f 2 K° 



G F 



E 



(5 pu a 2 + 2ajj + ^a^Ey, 
p 

r 3,EW " 

1 n ^,ti 1 



+ 



a 4 2 a 4,EW + ^3 2 ^3,EW 



(B24) 



(B25) 



V^B— = % E [W"l+A) + «S + <EW+^ + ^Ew]^ (B ' 

p=M,c t 

-V^o = % E A«[<5 pu /3 2 + ^-^< E w + /3 3 P + /33 P ew]^ ( ^ 2 ' F) ' 
^% a+ x- = ^ E A?[<5 P ^ + ^ + a^ EW + ^-i/3f )EW ]^ a2 ^ 

2 D=u . c z 



02,^) 



p=u,c 



B^— ►a !*' ~~ /o A P H "4 + 2 4,EW P3 + 2 ^ 3 ' EW J 

2 V ^ n=u.c I ~ 



p=u,c 



+ 



2 a 3,EW 



<BK,a 2 ) 



(B26) 
(B27) 
(B28) 



(B29) 



V2^ 



» p=u,c 



$pu «2 + 2«3 + -a| EW 



+ \/2 [<5 pu /3 2 + ag + - -a^ )EW - ^a£ EW + /f + /f jEW 



^(BX*,^'') 



«2 + 



C) 



+ 



<5pu (ai + ft) + + a^ EW + /3f + $ >EW X^.^) L 



\/2^ 



E a p s) < 



«2 + 2a§ + -Og EW 



+ ^[,5 pc a 2 + af]x( BX2 ^ 



+ 



"4 2 4 > EW ^ 3 2 3,EW 



Y( B vi'\K*2) 



(B30) 



(B31) 



and the amplitudes for i? — >■ If 2 71 " can be obtained from S — > Ka 2 with the replacement (K, a 2 ) 

(k;,7t). 
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2. B -^TV decays 



a. Decay amplitudes with AS = 0: 



The amplitudes for B — > fop can be obtained from B — > f2ir with the replacement (/2,7r) 

(f2, P ), 

Gi 



+ #* + - + ^ew]^ 1 '^ 

^ (a$ + ft) + 2af + of* + laf - iag w + + 2/f^ 



+ 



o^3,EW o^4,EW 



2 



p,h 
a 3,EW 



X (-B^/!)l (B32) 



2^4 



V ^ p=u lC L 

+ 2if B fiU 



rJ>> h - -rv p ' h 
"3 2 3,EW 



°4 2 4,EW 



°4 2 4,EW 



<fifi 



(B33) 



the amplitudes for B — > a^p can be obtained from B — > aiK with the replacement (02, 7r) — >■ (02, p), 



2 p= u ,c I L ^ 



p,h 
3,EW 



.In,?'' 1 , aP,h , oP,h 
2«4,EW ^ P3 + P 3> EW 



v (_Ba2,ai) 



+ 



5 pu (of + /3 2 ft ) + at + < w + + /f £ w 



■p,/i 1 ap^ 



X 



(Bui,d2) 



(B34) 



"2 = ^ E c ^ { [<W (<4 " tf) + 2a f + af h + lag £ w - ±aj* 



EW 



1 oP,h _ 1 op,7i ^ P^ 

" 2^ 3 ' EW 2 4 ' EW 



(Bo2,a;) 
ft 



+ 



<*pu (-«2 - + «4 ft - o^EW - o ff 



p,/i 



1 



~t~P3 2^ 3 ' EW 2 4 ' EW 



2"J,avv 2 ~~4,EW 



V ^ p= n ,c 



1 



p,h 



a 3 2 a 3' EW 



-V2^ a o, = ^ E ^ 

V p=U,C 



n 



1 

2 a 3,EW 



(B35) 
(B36) 
(B37) 



and the amplitudes for — >■ K 2 K* can be obtained from — >■ i^ 2 ^ with the replacement 
(K* 2 ,K) -»• (K2, K*) and c = -1. 
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b. Decay amplitudes with AS = 1 



The amplitudes for B — > fiK* can be obtained from B — V $2K with the replacement (/2, K) — > 
(f2,K*), and the ones for B — > K* 2 p and B — > K*a 2 can be obtained from B — > Ka^ with the 
replacement (K,a2) — > (K 2 ,p) and (K,a2) — > (K*,a2), respectively, 



y/2A* 



b-^k;- 



% E ^ 

V 2 p= u ,c 



3,EW 



(Bu;,K 2 ) 



h 



+ 



2 "3,EW 



-5° T7* 

B ->K w 



r h i r. p./i i 1 p,h 
<V«2 + 2«3 + -or' 

GF ^(s)} [_p,h 



X 



(BK 2 ,uj) 



~, E x p ] - 



O 



1 

2 a 4,EW 



+ 



SpuaZ + 2a p R ' h + 



1 



2 a 3,EW 



p,h 



x 



(bT 2 ,lo) 



p,h 



1 



p,h 



2 ^3,EW 



(Boj,K 2 ) 



(B38) 



(B39) 



.4 



.4^-0 — *0 



% E 4 S) + < ,h + h - > 

1 



1 p,h 

2 - J 3,EW 



--a 



p,h 
4,EW 



'3 P3.EW 



V * p=u,c 



c4' h + o% h 



" 2 a 3,EW ~~ 2 a4 ' EW 



"TP3 2^ 3 ' EW 



(B40) 



(B41) 



Appendix C: Explicit expressions of annihilation amplitudes 



The general expressions of the helicity-dependent annihilation amplitudes are given in Eqs. 
(f90|) - ff99]) . They can be further simplified by considering the asymptotic distribution amplitudes 
for <&y, & v , $ T and 3> t : 

$l ± (u) = 6uu, $„(«) = 3 (2u- 1), 

$jT ± (u) = 30uu(2u - 1), $ t (it) = 5(1 - 6u + 6u 2 ), 



We find 
,/,o 



4' -(VT) 



l $f (u) 



Jo 



-•'JO \/ -7ra s 
10 



6Xf - 23X° A + 22) + -(2X° A - 1)(X° A - 3)r 



^2 



vra. s 



(ex 



/ - 23X7, + 17) + 9(2X A - 3)(X A - 2)^1 



ray 



4'° (TV) = -A f 3 <°(VT), 4'-(TV) = -4'-(yT) 1 



3U\/ g 71 " 1 ^ 



30 



7ra, 



"3 (X A 2 - AX A 



X A ~ 2X A 



TUT 



7T 



2X\-Q + -\r x 



2)^ + (3X; 

my x V 



12X7 + 2tt' 



2 ^ my t 
ttit x 



(CI) 

(C2) 

(C3) 
(C4) 
(C5) 

, (C6) 
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30\ -ira s 
3 



4°(T7) = A%°(VT), 
A\ G {VT) 

A\ (TV) 

A%°(VT) = -A^°(TV), 
for TV modes, and 



A? 
3A^ + 4 



(TV) 

2 



A 1 ' 

^3 



(VT), 



^ + t(X\-2,)(X\-2) 



V T 

r x r x 



-30 



X A + 29-3^ + ^X A -3)(X A -2)r^ 



A^°(TV) 



-A\\VT), 



4i.PT) 
Ai.PT) 

A\(PT) 

A\(TP) 

4(PT) 
Al(TP) 



ira. 



WX A (6X A - ll)rj + 45(2X A - 1)(X A - 3)r 



na. 



30 



10a/ — vra, 
3 



AV - \X A + 4 + 



7T 



2X A - 6 + 



7T 



3(3A A + 4 - 7T 2 ) + -A A (X A - 3)r£r£ 



-10 

4(TP), 
4(iT), 



ira. 



3{X A + 29 - 3vr 2 ) + ^X A (X A - 3)r£r£ 



A£(TP) = A£(PT), 



4(PT) = 4(TP), 



(C7) 
(C8) 

(C9) 
(CIO) 

(Cll) 
,(C12) 

(C13) 

(C14) 

(C15) 
(C16) 



for TP modes. As pointed out in 63], since the annihilation contributions A\ 2 are suppressed 
by a factor of mir^/mfj relative to other terms, in numerical analysis we will consider only the 
annihilation contributions due to .A3' , A3' - , A^'° 3 and Ag' - . 

The logarithmic divergences that occurred in weak annihilation in above equations are described 
by the variable X A 



du 
u 



X A , 



1 in 77 
du 

u 



-AXa) 



(C17) 



Following 241 ] . these variables are parameterized in Eq. (llOOp in terms of the unknown real pa- 
rameters p A and 4>A- For simplicity, we shall assume in practical calculations that X\ are helicity 
independent, X A = X\ 



X\. 



Appendix D: The 77 — rj 1 system 

For the 77 and 77' particles, it is more convenient to work with the flavor states qq = (uu+dd)/\/2, 
ss and cc labeled by the rj q , rj s and rj®, respectively. Neglecting the small mixing with 77°, we write 



1 77) \ _ / cos0 -sin^A / |t7 9 ) 
1 77') / V sin<^> cos0 J \ \rj s ) 



(Dl) 



where <j) = (39.3 ± 1.0)° [65[ is the 77 — 77' mixing angle in the 77^ and rj s flavor basis. Decay 
^C')> and %o 



constants f q ( , v / s ( /) and f c (l) are defined by 



(Ol^/iT*?!^) = i-^f^Pv, (0|s7/.75sh (,) ) = if*(>)Pn, (0|c7a»76c|t/^) = if°(>)P», (D2) 



6 A different mixing angle = (35.9 ± 3.4)° was obtained recently in the analysis of 66] based on vector 
meson radiative decays. 
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while the widely studied decay constants f q and f s are defined as [651 ] 



(0|s7^7 5 s|r? s ) = ifsP 11 . 



(D3) 



The ansatz made by Feldmann, Kroll and Stech (FKS) 65J is that the decay constants in the quark 
flavor basis follow the same pattern of 77 — 77' mixing given in Eq. (|D1|) 



fq fs 

J Tj Jr\ 

fl fs 

J„' in 1 



COS <p — SUl ( 

sin 6 cos d> 



fq 
fs 



(D4) 



' rj' ■> r\' 

Empirically, this ansatz works very well [651 ]. Theoretically, it has been shown recently that this 
assumption can be justified in the large- N c approach [671 ]. 

It is useful to consider the matrix elements of pseudoscalar densities 681 ] 

2m q (0\qj 5 q\rj^'' - — 

and define the parameters h q and h s in analogue to f q and /. 



/gfcjoi 2m s (0|s7 5 s|?7 ( ' ) ) = ih 8 ^, 



(D5) 



1 

2m q (0\qj 5 q\r]g) = -j=h q , 2m s (0|s7 5 s|?7 s ) = ih s , 



(D6) 



and relate them to by the similar FKS ansatz as in Eq. ([D4|) 



In this work, we shall follow 



571 ] to use 



h> = 0.0013 GeV 3 , 



-0.0555 GeV 3 , = 0.0011 GeV 3 , 



/« = 109 MeV, 



-2.3 MeV, 



•111 MeV, 
-5.8 MeV, 



/J = 89 MeV, 



h s rf = 0.068 GeV J 
= 136 MeV, 



in 



'hi 



741 MeV, 



in 



where we have used the perturbative result 69] 



(') 



I2m 2 c y/2 



802 MeV, (D7) 



(D8) 



The form factors for B — > rj^ transitions obtained in QCD sum rules are 70j] 



F ^(0) = 0.229 ±0.024 ±0.011, 
F^'(0) = 0.188 ±0.002t5|± 0.019 ±0.009, 



(D9) 



where the flavor-singlet contribution to the B — > rfi^ form factors is characterized by the parameter 
Bn, a gluonic Gegenbauer moment. It appears that the singlet contribution to the form factor is 



small unless -Bf assumes extreme values ~ 40 [701 ]. Using the relation 



Fo,i = F oT cos & F o7 = F oT sm h (D10) 

we obtain Fq^ 9 (0) = 0.296 ±0.028 as shown in Table ITVl The momentum dependence of the form 
factor can be found in [7o| . 



jBlJg 
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